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An alternative approach - nonequilibrium evolution thermodynamics, is compared with classical 
Landau approach. A statistical justification of the approach is carried out with help of probability 
distribution function on an example of a solid with vacancies. Lwo kinds of kinetic equations are 
deduced in terms of the internal energy and the modified free energy. 
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A great number of works devoted to development of 
the phase-field approaches for description of different na- 
ture phenomena has been published recently [H-Q . They 
are continuing the Landau theory of phase transitions 
[Tol [TT| . During the last decade an alternative approach, 
close to the Landau technique of phase transitions, has 
been developing too [12h15| |. There is no any statistical 
justification of the approach and in the present article 
such justification is carried out with help of probability 
distribution function (PDF) on the example of a solid 
with vacancies. In the approach the conjugated pair 
of thermodynamic potentials, that is, the internal and 
(modified) free energy, is introduced. The main intrigue 
is that the equilibrium or stationary state does not co- 
incide with extremes of such thermodynamic potentials 

mill- 

If a solid consisting of N particles has n point defects 
(e.g., vacancies, substituted atoms, etc.) then in this 
case the equilibrium (or stationary) state can be found 
from the maximum of PDF taken in the form [n], Il6l - fl8j 
(curve 1 in Fig. [1]) 



f(n) = CWexp(- 



U(n) 
kT 



(1) 



where C is a normalizing factor, U(n) is the internal en- 
ergy dependent on the number of structural defects, k is 
the Boltzmann's constant and W is the thermodynamic 
probability as a possible combination of the number of 
defects n and the number of atoms N of a crystal [l|| 



W = 



(N + n)l 
N\nl ' 



(2) 



the logarithm of which is a configurational entropy (curve 
2 in Fig. 1) 



S c = khxW. 



(3) 



Note, that the configurational entropy is a one-valued 
function of the number of defects. It is absolutely in- 
dependent of the energy of defects (and of temperature 
too). 




n(%) 



FIG. 1: Basic set of thermodynamic functions for quadratic 
approximation: 1 - probability distribution function /; 2 - 
classic configurational free energy F c ; 3 - internal energy U; 
4 - a tangent to U at the equilibrium point; 5 - modified 
configurational free energy F c ; 6 - effective internal energy U; 
7 - effective configurational free energy F c . The right-hand 
scale is for PDF. 



The pre-exponential factor of (JTJ describes the combi- 
national, that is entropic, part of the distribution func- 
tion, related to degeneration of macrostates. The expo- 
nential factor describes the restrictive part of the distri- 
bution function, connected with the overcoming of po- 
tential barriers between microstates. In a power presen- 
tation 



K 



U = U +J2 



(-1) 



k-l 



-Uk-in 



(4) 



fc=i 
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where U is the internal energy of a solid with defects, 
Uq is the internal energy of the defect- free solid, Ufe_i 
are some coefficients, K is the highest power taken into 
account. In a linear approach U — Uq + UqU the same 
potential was introduced by Frenkel [16j . In a quadratic 
approximation the graph of U is shown in Fig. [TJ (curve 
3). 

Bringing the variables independent of the number of 
defects n into the "inessential" constant C we can write 



2 



the expression ([lj in the form 



f{n) = C : exp( — — ), 



or as a product 



JV 



/(») = C JJ(n + i) exp(- 



kT 



U{n) 
~k~T~ 



Through differentiation we obtain 



df(n) 
On 



N-\-n _^ n 1 



fc=l 



where 



= d[/(n) 
dn 



fc=i 



fcT 



)/(n), 



^(-i)W 



(5) 



(6) 



(7) 



(8) 



fe=0 



is the energy of a defect (tangent 4 to the internal energy 
in Fig. [TJ. 

The extreme value of probability distribution function 
is at n, it obeys the following transcendental equation 



7V+n ^ n 1 



fc=l 



fe=l 



0. 



(9) 



The first two terms are sums S of the a slowly divergent 
harmonic series. This part of the equation depends only 
on system size and doesn't depend on material parame- 
ters. It is the fundamental part of Eq. [5] decreasing with 
the growing parameter n (curve 1 of Fig. [2]). The last 
term in Eq. |9] depends on material parameters through 
coefficients Uk . At different relations between these coef- 
ficients Eq. [S] may have several solutions. For example, 
in cubic approximation, Eq. [9] may have sole solution 
(the uppermost curve and the two lower curves of series 
2 in Fig. [5J or two stable solutions (the two intermediate 
curves) with dependence on parameter uq. To calculate 
the fundamental curve S the value N = 2000 was chosen. 
As seen, the positions of routes of Eq. [9] coincide with 
maxima of PDF calculated directly by Eqs. (d])-®. 

The second upper curve of series 2 in Fig. [2] corre- 
sponds to coefficient uq of numerical value ~ 0.15eV". For 
many real materials the energy of vacancy formation is 
higher, for example, in the case of cooper it is about leV, 
where Eq. [9] has the single solution. Nevertheless, this 
coefficient can be reduced at the cost of tensile stresses 
[TJ, [20j], size effects [HI or a temperature rise [H^, 
and, as a result, a region with two solutions of Eq. [S]can 
be reached. In this case, a structural phase transition 
between two stable states, one of which is in the region 
of low defectiveness, and the other in the region of high 
defectiveness, is possible. The kinetics of such structural 
phase transition will be considered below. 

From the table value of partial sums one can find (25j | 



1 

k=l 



In n 



1 

2n' 



(10) 




FIG. 2: Dependence of fundamental curve S (1), en- 
ergy of vacancy u (2) and PDF (3) on the number 
of vacancies n at different values of parameter uq : 
0.135; 0.13; 0.125; 0.12; 0.125(e , l/). The rest parameters are 
mi = 0.6 • W~ s (eV), u 2 = 0.125 ■ 10 _5 (eV), T = 300^. The 
right-hand scale is for PDF. 



where C is a constant. By substituting (fTO)) into ([9]), 
for TV >> n >> 1 we obtain a relation between the 
number of defects n and average energy of defect u in 
the equilibrium state 



n = N exp(- 



(11) 



As evident from the last formula, the energy of defect is 
not strictly constant, but it depends on the total number 
of defects. Relation (jlll) is an equation of state for the 
equilibrium case, relation ([8]) is the equation of state too, 
but for a more general nonequilibrium case including the 
equilibrium state as a partial case. Eqs. ([8]) and (|TT|) need 
be considered together as a set of equations for deducing 
both the energy of defect u e and the density of defects 
n e in the equilibrium state. 

Thus, the equation of state (fTTj) is obtained from con- 
dition of the most probable state as the maximum of 
probability distribution function ([T]). Same result can be 
obtained from the principle of free energy minimum. Re- 
ally, the above-described procedure can be schematically 
displayed as [TjJ 



U(n) - kTXnW = U ~TS C = F c 



(12) 



It need be mentioned that product TS C introduced into 
the definition of the canonical free energy F c is a bound 
energy lost by the system for doing work. On the other 
hand, the total energy of defects in the main part is physi- 
cally the energy lost for the production of work too. Only 
a little part of it remains for the work production. It fol- 
lows that 



TSr. 



(13) 



And now we can introduce a new specific kind of the 
free energy by subtracting the bound energy in the form 



3 



of product un from the internal energy ((4]) . In quadratic 
approximation (Fig. [1]) 



the "chemical potential" is in the form 



F c = U - un = U + —{u - u) 2 . 

2u\ 



(14) 



Here we use the equation of state © to eliminate the 
density of defects. It is easily found that 



dF c 
du 



(15) 



The both relations flSJ) and (IT5l) are a connected cou- 
ple of equations between the internal energy U and the 
modified configurational free energy F c on one hand, and 
between the density of defects n and the energy of defect 
u on the other hand. One can see that the density of 
defects is the eigen-argument for the internal energy, and 
the energy of defect is the eigen-argument for the mod- 
ified configurational free energy. Hence, the exact free 
energy F c is expressed according to © and (|T2"j) through 
variable n, which is not its eigen-argument. 

Because the energy needed for the formation of a new 
defect is smaller, in the presence of others than, in defect- 
free crystal, the quadratic term in Q has negative sign. 
Note that expression ^ is true for both equilibrium and 
non-equilibrium states. In this approximation the inter- 
nal energy is a convex function of the defect number hav- 
ing the maximum at point n — n max , as shown in Fig. 
[T] a. In the same approximation the modified configura- 
tional free energy is a concave function with the minimum 
at point u = uq. 

With relationships ([5} and ([TT| it is easy to show that 
the stationary state corresponds neither to the maximum 
of the internal energy nor to the minimum of the free 
energy F c . The stationary state is at point n = n e , where 



dU 
dn e 



dF c 
du e 



(16) 



Here the additional subscript e denotes the equilibrium 
value of a variable. 

If the system has deviated from the stationary state it 
should tend back to that state at a speed the higher, the 
larger the deviation [l3l - fl5| 



dn 



T7«( 



dF c 



(17) 



du 

dt ,ay du 

The form of kinetic equations (|17l) is symmetric with 
respect to the use of internal and configurational free en- 
ergy. If the equilibrium state is closer to the maximum 
of internal energy, then in view of stability of the solu- 
tion the upper sign is chosen (convex function), if it is 
near the minimum of internal energy, the lower sign is 
taken (concave function). The both variants of the ki- 
netic equations are equivalent and their application is a 
matter of convenience. 

In the right side of the well-known Landau- 
Khalatnikov kinetic equation [l(| 

dn dF c 



dt 



-r 



dn 



(18) 



dF c 
dn 



(19) 



From the thermodynamic point of view, a variable of 
such kind is not chemical potential really, as is specified 
by the "argument" foreign to the free energy. Still, this 
notion can be used in practical work, as it directly realizes 
the minimization principle for the canonical free energy. 

If we assume that the equilibrium energy of defect u e 
and the number of defects n e are slowly changing during 
an external action then their product can be introduced 
under differentiation sign in (|17[) and a new kind (shifted) 
of internal and free (effective) energies can be defined 



U = U — u e n, F c = F c + un e . 
Then equations (|17[) are simplified a little 



(20) 



dn dU du dF c 

m =±ln ^ at = ^^7" (21) 

The original potentials U and F c are connected by 
means of a Legendre-like transformation 



F c = U — un. 



(22) 



The shifted potentials U and F c are connected by 
means of transformation 



U — un + u e n + un e , 



(23) 



which differs from the Legendre-like transformation by 
anticommutation bracket [un] = u e n + un e . 

Graphs of the original and effective internal energies 
are given in Fig. [3] In the considered interval the orig- 
inal internal energy has no extremes, but its curvature 
changes the sign. In the region of a small number of 
vacancies it is convex, for a larger number it is concave 
(curve 2). That's why the effective internal energy has 
maximum in the first case and minimum in the second. 

The stationary point for the shifted potentials coin- 
cides with maximum of U (Fig. [3]) and with minimum 
of F c for the upper signs in (|2Tj) . Thus, U is an effective 
thermodynamic potential, for which the tendency of the 
original part of internal energy to minimum is completely 
compensated by the entropic factor. Twice modified con- 
figurational free energy F c tends to minimum, but this 
tendency differs from the case of the canonical config- 
urational free energy F c . The effective thermodynamic 
potential F c tends to minimum in the space of eigen- 
argument u, while the canonical free energy F c tends to 
minimum in the space of non-eigen "argument" n. 

In the paper, a phenomenological approach based on 
generalization of the Landau technique is considered. For 
fast processes thermal fluctuations have no time to exert 
essential influence and it becomes possible to consider 
the problem in the mean-field approximation. The ap- 
proach is based not on an abstract order parameter but 
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on physical parameters of structural defects - their quan- 
tity (density) and average energy. The new more gen- 
eral form of kinetic equations, symmetric with respect 
to using the internal energy U and the modified config- 
urational free energy F c , is proposed. In this case, the 
density of defects and defect energy are related by sym- 
metric differential dependences of type (jSJ , (fT5|) and (|16[) . 
As the defect energy in the stationary state is not zero, 
the extreme principle of equality to zero of the deriva- 
tive of free energy with respect to "order parameter" 
in the framework of nonequilibrium evolution thermody- 
namics breaks down. This principle need be substituted 
with the principle of the tendency to a stationary state. 
Stationary-state characteristics can not be determined in 
the framework of phcnomcnological approach, statistical 
and microscopic approaches are required. 

The present form of kinetic equations can be general- 
ized to all types of regular or randomly distributed de- 
fects. 
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FIG. 3: Original (a) and effective (b) internal energies as 
functions of the number of defects: a) 1 is family of original 
internal energy; 2 is linear transform of internal energy with 
same curvature; b) 1 — 5 are positions of the equilibrium states 
for the same parameters as in Fig. [2] 



The work was supported by the budget topic 
0106U006931 of NAS of Ukraine and partially by the 
Ukrainian state fund of fundamental researches (grant 
F28. 7/060). 



[1] I. S. Aranson, V. A. Kalatsky, and V. M. Vinokur, Phys. 

Rev. Let. 85, 118 (2000). 
[2] A. Karma, D. A. Kessler, and H. Levine, Phys. Rev. Let. 

87, 045501 (2001). 
[3] L. O. Eastgate, J. P. Sethna, M. Rauscher, T. Cretegny, 

C. S. Chen, and C. R. Myers, Phys. Rev. E 65, 036117 

(2002). 

[4] V. I. Levitas, D. L. Preston, and D. W. Lee, Phys. Rev. 

B 68, 134201 (2003). 
[5] J. C. Ramirez, C. Beckermann, A. Karma, and H. J. 

Diepers, Phys. Rev. E 69, 051607 (2004). 
[6] L. Granasy, T. Pusztai, G. Tegze, J. A. Warren, and 

J. F. Douglas, Phys. Rev. E 72, 011605 (2005). 
[7] C. V. Achim, M. Karttunen, K. R. Elder, E. Granato, 

T. Ala-Nissila, and S. C. Ying, Phys. Rev. E 74, 021104 

(2006). 

[8] A. Svandal, Modeling hydrate phase transitions using 
mean-field approaches (University of Bergen, Bergen, 
2006). 

[9] J. Rosam, P. K. Jimack, and A. M. Mullis, Phys. Rev. E 
79, 030601(R) (2009). 
[10] L. Landau and I. Khalatnikov, Sov. Physics Doklady 96, 
459 (1954). 

[11] L. Landau and E. Lifshitz, Statistical mechanics (Perga- 

mon Press, Oxford, 1969). 
[12] L. S. Metlov, and A. F. Morozov, High Pressure Physics 

and Technics 7, 3, 58 (1997). 



[13] L. S. Metlov, Bulletin of Donetsk Univ. Ser. A: Nature 

Sciences, 1, 169 (2007). 
[14] L. Metlov, Bulletin RAS, Physics 72, 1283 (2008). 
[15] L. Metlov, Phys. Rev. E 81, 051121 (2010). 
[16] Ya. I. Frenkel, Kinetic Theory of Liquids (Dover, New 

York, 1955). 

[17] A. S Steinberg, Reportage from alloys world (Phys- 
MatLit, Moscow, 1989). 

[18] Yu. M. Gufan, Soros education journal 7, 109 (1997). 

[19] R. W. Balluffi, R. W. Siegel, K. H. Lie and D. N. Seid- 
man, Determination of concentration and formation en- 
ergies and entropies of vacancy defects from quenching 
experiments. Report No 125 (State university of New 
York at stony brook, New York, 1969). 

[20] S. Chaudhry and M. E. Law, J. Appl. Phys. 82, 1138 
(1997). 

[21] W. H. Qi and M. P. Wang, Physica B 334, 432 (2003). 

[22] J. M. Richan, R. Najafabadi, L. Zhao, and 
D. J. Srolovitz, J. Phys.: Condens. Matter 1, 539 (1992). 

[23] L. Zhao, R. Najafabadi, and D. J. Srolovitz, Modelling 
Simul. Mater. Sci. Erg. 4, 4923 (1993). 

[24] A. I. Karasevskii, and V. V. Lubashenko, Phys. Rev. B 
71, 012107 (2005). 

[25] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Se- 
ries, and Products ( Academic Press, New York, 2007). 



